Dynamical-charge neutrality at a crystal surface 
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For both molecules and periodic solids, the ionic dynamical charge tensors which govern the infrared 
activity are known to obey a dynamical neutrality condition. This condition enforces their sum to 
vanish (over the whole finite system, or over the crystal cell, respectively). We extend this sum rule 
to the non trivial case of the surface of a semiinfinite solid and show that, in the case of a polar 
surface of an insulator, the surface ions cannot have the same dynamical charges as in the bulk. The 
sum rule is demonstrated through calculations for the Si-terminated SiC(OOl) surface. 
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The basic quantity addressed in this work is the dy- 
namical charge of a given ion s in different environments: 
within a molecule, in the bulk of a crystalline solid, and 
at a solid surface. This charge is a cartesian tensor, Z*, 
which has the point symmetry of the ionic site: its com- 
ponents Z* measure the dipole linearly induced (in 
the a direction) by a unit displacement of the ion s (in 
the (3 direction) Equivalently, Z* measures the force 
linearly induced on the given ion s by a unit electric 
field (at zero displacement): the dynamical charges gov- 
ern therefore the infrared activity of the system. The 
two cases of molecules [D and of bulk solids [^jH have 
received previous attention in the literature, both as a 
matter of principle and as a subject of practical calcula- 
tions, while the case of a crystal surface has never been 
considered In a neutral molecule the sum (over all 
the ions) of the dynamical charges must vanish, since 
a rigid translation of the molecule as a whole induces 
no dipole: we will refer to this sum rule as to dynam- 
ical neutrality. The analogue in a crystalline dielectric 
goes under the name of acoustic sum rule (ASR), and is 
spelled out in an equally simple manner: the dynamical 
charges sum to zero over the crystal cell. However, the 
underlying theory is by far less trivial ||^. We show here 
that a crystalline surface must be dynamically neutral in 
order to ensure that a rigid translation of the semiinfinite 
solid as a whole does not affect the work function. Ex- 
plicit formulation of such neutrality requires in general 
the regularization of a nonconvergent sum. The result is 
a constraint for the dynamical charges of the surface ions: 
in the particular case of a polar surface — such as (001) in 
the zincblende structure — the novel sum rule forbids ions 
of a given chemical species to have the same dynamical 
charge at the surface and in the bulk. 

The dynamical charges — both in molecules and in 
solids — have in general nothing to do with the nomi- 
nal static charge of the given ion. The difference be- 
tween dynamical and static charges is particularly dra- 
matic when the material has a mixed ionic-covalent char- 
acter [P: essentially, the dynamical charge measures the 



current flowing along the bonds when the bond-lengths 
are varied. In the molecular case the charge perturba- 
tion induced by an ionic displacement is localized, and 
the deflnition of Z* is straightforward and unique 
In a crystalline solid one refers instead to the rigid dis- 
placement of a whole sublattice: since the induced charge 
is a periodic function, one has to carefully specify the 
boundary condition assumed in the solution of Poisson 
equation. Usually one defines two kinds of tensors: the 
transverse (or Born) charge 2.*s^'^\ and the longitudinal 
(or Callen) one Z^^'"''. The former is the most fundamen- 
tal when deahng with bulk properties jsj: it is defined 
via the macroscopic polarization AP linearly induced in 
the solid by a rigid displacement of the s sublattice by 
an amount u^. , while the field is kept vanishing: 

(1) 

where 51 is the cell volume. The longitudinal charge Zs 
is analogously defined, but the sublattice displacement is 
performed in a depolarizing field AE = — 47rAP. The 
relationship between the two is Z^*''^^ = Soo'Z-l^^\ where 
Eoo is the electronic dielectric constant (we assume it 
isotropic for the sake of simplicity). The sum rule is 
usually stated for the transverse charges |^, but equiv- 
alently holds for the longitudinal charges as well. We 
wish to deal with molecules, surfaces, and bulk solids all 
on the same footing: then the longitudinal charge is the 
quantity of choice in the crystalline case. Suppose in fact 
we displace by an amount only a single ion of species s 
in an infinite — and otherwise unperturbed — crystal: the 
induced charge is localized [Q, and it is easy to prove 
that its dipole d is (to linear order) 

rf"=ECaK/^- (2) 

This can be regarded as an equivalent definition of the 
longitudinal charge. An identical definition as Eq. (||) 
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holds for the (unique) dynamical charge in the molecular 
case. We therefore drop the superscript (L) in the follow- 
ing, and we assume Eq. (||) as a uniform definition of the 
dynamical charge for a molecule, a bulk solid, and a sur- 
face: whenever a bulk charge is concerned, wc implicitly 
refer to the longitudinal one. 

We are now ready to discuss the surface problem. Sup- 
pose we have a semiinfinite crystalline insulator, where 
we distinguish a surface region and a bulk region. We 
assume that the macroscopic field vanishes both in the 
bulk of the solid, and outside in the vacuum region: this 
amounts to require that no static charge is present at 
the surface. In these hypotheses, the average of the elec- 
trostatic potential in the bulk of the solid with respect 
to the vacuum level is a well defined quantity, which in 
fact determines the work function of the given surface. 
We are going to impose the physical requirement that 
the work function is not affected by a rigid translation of 
the semiinfinite solid as a whole: a necessary condition 
for this to occur is our novel sum rule for the dynamical 
charges. 

Suppose the surface is normal to the z axis: the system 
has then a two-dimensional periodicity normally to z: we 
indicate with A the area of the unit cell. We consider 
a rigid displacement of an ionic layer, i.e. we displace 
by an amount Ug a given ion and all its translationally 
equivalent ones: this rigid displacement induces a dipole 
per unit area, hence a potential lineup across the layer 
whose value is given by: 

^'i>=^^^K.^^i■^s,^}■ (3) 

A rigid translation of the semiinfinite crystal by an 
amount u induces a total lineup which is, by linearity, 
the sum over s of the expressions in Eq. (|^) . As antici- 
pated above, we explicitly require this lineup to vanish, 
hence a naif expression for the constraint appears to be: 

5] Z* 3^ = (any/?). (4) 

s 

This formal expression for the dynamical charge neutral- 
ity of the surface cannot be used as such, given that the 
infinite sum in general does not converge: in the bulk 
region it oscillates periodically, owing to bulk dynamical- 
charge neutrality. 

Therefore the problem is to regularize the indetermi- 
nate sum in Eq. by using the appropriate physical cri- 
terion. The universal panacea for this class of problems 
at large is the macroscopic average introduced in Ref. ||] 
and widely used by several authors: by construction, 
the macroscopic average yields the correct electrostatic- 
potential average in the bulk region. Application to the 
present case is straightforward. One first maps the prob- 
lem into a simple electrostatic one by assigning a point 
charge of magnitude Z*3^ to each ion. Secondly, one 



evaluates the planar average of this charge distribution, 
which takes the general form: 

p(z)^i^Z;3^5(z-z,), (5) 

s 

where Zs are the positions of the ionic planes along the z 
axis. Finally, one filters p(z) through the convolution: 

1 j-z+b/2 

=p{z) = T dz'piz'), (6) 

Jz-b/2 

where b is the one-dimensional periodicity of the bulk 
region. The result is a piecewise constant function, van- 
ishing both in the vacuum and in the bulk, and whose 
integral in the surface region is trivial. We state the sum 
rule by requiring this integral to vanish. It is important 
to realize that — at variance with the naif expression of 
Eq. (^) — the surface sum rule involves both the values of 
the dynamical charges and the coordinates of the planes 

Zs- 

The nontrivial content of the sum rule will be made 
clear with two examples. We consider simple polar sur- 
faces of a cubic binary crystal, whose bulk dynamical 
charges have the form Z*^^^""^ = i-iy\Z*\ S^t}. The 
/? ^ 3 components of the sum rule expressed by Eq. (^) 
are correct, because of periodicity, for the same reason 
why this equation is correct in the bulk. We focus then 
on the (3 = 3 component, and we further limit ourselves 
to the cases where the ionic planes parallel to z contain ei- 
ther cations only, or anions only, alternately. This is the 
case for the (001) and (111) surfaces in the zincblende 
structure, which we are going to illustrate separately, as- 
suming an ideal (truncated-bulk) geometry i.e. neglect- 
ing surface relaxation and recostruction. 

We start with the (001) surface, where b is one half of 
the cubic lattice constant and the ionic planes are equally 
spaced by an amount 6/2. We assume that the crystal 
lies in the positive-z half-space, so that the coordinate 
of the s-th plane can be assumed to be (s — 1)6/2. The 
piecewise constant function 'p{z) vanishes from — oo to 
—5/2; it assumes the value .^133 (apart for a constant 
factor) up to z = 0; then ^(z) = ZJ' 33 -I- Z2 33 up to 6/2, 
and in general Z* + Z*^^ 33 in the following intervals, 
until it vanishes again in the bulk region, say for s > N. 
Since the the intervals are of equal length, the sum rule 
for the integral of p(z) requires: 

N 

2 ^ -^3,33 + -^^+1,33 = 0- (''') 

We arrive thus at an outstanding finding: the dynamical 
charges of the ions in the surface region sum up to one 
half of the bulk dynamical charge (with the appropriate 
sign). This is precisely the relationship that one would 
obtain by replacing the dynamical charges with static 
point-like charges and requiring the surface to be neutral. 



2 



Switching now to the (111) surface, the linear period 
b is \/3 times the lattice constant, and the ionic planes 
are no longer equally spaced: the sum of the dynamical 
charges in the surface region is then a certain fraction 
(±1/4 or ±3/4) of the bulk dynamical charge. For any 
polar surface we have therefore a strong constraint, which 
e.g. forbids the surface ions to have the same dynamical 
charges as in the bulk. 

We illustrate our main finding using first-principles cal- 
culations of the surface dynamical charges for a paradig- 
matic test case: to keep matters simple we deal with an 
insulating surface. Our choice is the (001) Si-terminated 
surface of the zincblende semiconductor SiC, which has 
been the subject of recent theoretical work |]9|,|l0|. The 
actual structure is 2x1 reconstructed, but even the ideal 
(truncated bulk) one is insulating, and fits well our pur- 
pose of dealing with a test case as simple as possible: the 
sum rule for this case takes precisely the form of Eq. (0). 
All calculations are performed using density-functional 
theory in the local-density approximation ||l^. Most 
technical ingredients are pretty standard: plane-wave 
basis sets and norm-conserving pseudopotentials [ p^ , 
Ceperley- Alder exchange-correlation potential ||l^, and 
special-point Brillouin-zone sampling. We use plane 
waves up to a kinetic energy cutoff of 18 Ry. This ba- 
sis set is too small if one aims at a precise prediction 
of the physical properties of any carbon-based materials, 
but it is enough to demonstrate the points of principle 
addressed here, because the ASR in the bulk is well sat- 
isfied. We use a set of 28 irreducible special points in the 
bulk calculations, and a consistent set in the supercell 
ones. This large number of special points is needed to 
satisfy the ASR |5|. 

Our calculations provide a theoretical lattice constant 
of 8.24 a.u. (Expt. 8.25). Density-functional perturba- 
tion theory (DFPT) [|5|, implemented as in Ref. [p^ , 
provides for the bulk solid Zg^ = +2.58, Z^^^^ ^ 
—2.57, and Soo — 7.40. The experimental estimate [Q is 
|^*(T)| _ 2.7. Notice that the algorithm does not enforce 
neutrality: the error monitors the numerical accuracy of 
the calculation. The longitudinal dynamical charge of 
our bulk material is therefore \Z* \ = 0.34. 

The surface calculations have been performed in a su- 
percell geometry, where the SiC slab has a double Si ter- 
mination. We plot in Fig. |l| the pseudopotential valence- 
electron density for the 9-atom supercell; we also plot 
the macroscopic average, Eq. (|[), of the same density. 
Fig. clearly shows that the two surfaces are well sepa- 
rated, and that the central region of the slab is bulklike, 
with eigth electrons per cell in average. We have also 
accurately checked that each of the two Si-terminated 
surfaces is statically neutral. 

Starting from the reference equilibrium ground state 
of Fig. |l|, we calculate the dynamical charges of all the 
ions in the supercell using again DFPT for our composite 
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FIG. 1. Self-consistent valence electron density in the 
9-atom computational supercell, in units of electrons per bulk 
cell. Vertical bars on the abscissae indicate the positions Zs of 
the ionic planes. The solid line shows the average over planes 
parallel to the surface, as a function of the normal coordi- 
nate z. The dashed line shows the same function after the 
macroscopic-average filtering. 

structure of periodically repeated slabs. Obviously, the 
calculated macroscopic dielectric tensor is a rather ar- 
tificial quantity, only indirectly related to genuine mate- 
rial properties; equally artificial are the transverse charge 
tensors. As already stressed, the longitudinal dynami- 
cal charges are the relevant physical quantities in this 
problem: they are anisotropic tensors in the surface re- 
gion, and they must converge to their (isotropic) value 
±0.34(5q,/5 in the bulk region. The convergence proves to 
be rather slow: insofar as the dynamical charges are con- 
cerned, the surface region is much larger than the equi- 
librium charge density of Fig. || would suggest. Despite 
this fact, our slab is large enough to recover the bulk Z33 
value for the center Si ion (0.33 vs. 0.34 from the bulk 
calculation). The calculated relevant effective charges 
are: Zl,,{St) = +0.13, Zl,,{C) - -0.27, ^3*33(5*) = 
+0.29, 2'|,33(C) = -0.31, 2'5,33(5i) = +0.33. The ro- 
bustness of these figures has been then checked in two 
different ways. Firstly we have performed similar calcu- 
lations on a fully relaxed (though unreconstructed) struc- 
ture: the outermost SiC bond length increases by 4%. 
The picture is unchanged, but the bulklike limit is re- 
covered faster, with Z33 = 0.34 for the center Si ion. 
Secondly the adequacy of our 9-atom supercell has been 
checked against a few test calculations performed with 
a slab of 13 atoms. Using the data reported above, the 
sum of the dynamical charges Z* 33 over five layers in the 
surface region, up to the central Si, is 0.17, thus demon- 
strating our main finding of Eq. (|^). 

A moment of reflection shows that the surface sum rule 
looks like equivalent to the previously known (molecular) 
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sum rule for the finite slab. However, this is a very special 
case, where two equivalent polar surfaces are possible. In 
a more general case — such as SiC(lll) — terminating the 
bulk with two equivalent surfaces is impossible. Wc have 
then a different sum rule for each surface of the slab: only 
the sum of the two can be related to a previously known 
sum rule. 

In this work we have discussed the concept of ionic dy- 
namical charge at the surface of a crystalline dielectric. 
Wc have shown that the longitudinal charge is the most 
fundamental quantity, at variance with bulk problems 
where the transverse one plays the major role. A neutral 
surface must also be dynamically neutral, and this im- 
poses a nontrivial constraint for the dynamical charges 
near a polar surface. We have performed a case-study 
calculation for the Si-terminated SiC(OOl) surface. The 
results, besides demonstrating the novel sum rule, also 
show (in this particular case at least) that the dynamical 
charges converge to their bulk value more slowly than the 
equilibrium electronic density would suggest. 
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